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1. Introduction 


Nowacki, J., Panagiotopoulou-Stawnicka, H., Rutkowski, R. & Wasala, R. 
2019: Taxonomic separateness of the subspecies of Noctua interjecta Hiibner 
(Lepidoptera: Noctuidae) inhabiting central Europe. — Entomol. Fennica 30: 
33-42. https://doi.org/10.33338/ef.79904 


This study focuses on morphological and genetic differences between two cen- 
tral European subspecies of the noctuid moth Noctua interjecta Hiibner: N. 
interjecta interjecta Hübner, 1803 from Transcarpathia in Ukraine and N. inter- 
Jecta caliginosa (Schawerda, 1919) from Poland. While the morphological dif- 
ferences between these taxa are fairly obvious, individuals from these two popu- 
lations were found to differ also genetically, as indicated by the cytochrome c 
oxidase I (COI) gene sequences. The Transcarpathian population exhibits solely 
the NImizH3 haplotype, but the Polish one the relatively distant NImizH1 and 
NImizH2 haplotypes. A Neighbour-Joining tree and a network of haplotypes 
separated the two haplotypes found in Poland (N. i. caliginosa) from those de- 
rived from NImizH3 (N. i. interjecta). The divergence time for these two groups 
was estimated as 400—600 thousand years before present, coinciding with the 
maximum extent of the Pleistocene ice-sheet covering Europe. 
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southern part of its range, being distributed from 
Spain and Portugal across southern and central 


Noctua interjecta Hübner, 1803 is a moth species 
of the family Noctuidae with an Atlantic-Medi- 
terranean distribution. It currently occurs from 
the Iberian Peninsula across the whole of western 
and southern Europe reaching the Baltic Sea in 
the north (Fibiger 1993, Nowacki 1998, Aarvik et 
al. 2017). Only the nominative subspecies N. in- 
terjecta interjecta Hübner, 1803 is found in the 


France, south of the Alps across Italy and the Bal- 
kan Peninsula as far as Bulgaria. Noctua inter- 
Jecta interjecta is also found in central Europe 
south of the Sudetens and Carpathians, i.e. in 
Slovakia, Hungary, Romania and Transcarpathia 
in Ukraine. 

In contrast, the subspecies N. interjecta cali- 
ginosa (Schawerda, 1919) is found in the north- 


34 


ern part of the species’ range, from western 
France and the U. K. through Belgium, the Neth- 
erlands, Denmark, Norway and Germany to the 
countries around the Baltic Sea, i.e. southern 
Sweden and Finland, western and northern Po- 
land, and also Lithuania, Latvia and Estonia 
(Aarvik et al. 2017). 

Before the end of the 20" century, the only 
Polish records of N. interjecta were froma few lo- 
calities in the west of the country, from the Baltic 
coast as far south as Lower Silesia (Nowacki 
1994). These moths represented N. i. caliginosa. 
The last 30 years, however, have witnessed a 
rapid eastward expansion of the species within 
Poland and in other countries around the Baltic. It 
has been recorded in numerous localities from 
north-eastern Poland through central Wielko- 
polska (Greater Poland) to Upper Silesia (Blaik et 
al. 2009). In 2014, a single specimen of N. inter- 
Jecta was trapped for the first time in the Biesz- 
czady Mts. in south-eastern Poland (Nowacki & 
Patka 2015). Itis worth emphasising that, judging 
by its external appearance, this specimen un- 
doubtedly belonged to the subspecies N. i. inter- 
jecta. Hence, this is currently the northernmost 
locality of this southern form of N. interjecta in 
the eastern part of its central European range. 
Consequently, both subspecies of N. interjecta 
are found in Poland, albeit in localities separated 
from one another. 

This paper examines the hypothesis that the 
two central European subspecies of N. interjecta, 
at present occurring as geographically isolated 
populations and exhibiting marked morphologi- 
cal differences, are also genetically distinct. This 
may be the result of the originally contiguous 
range in Europe having become disjunct in the 
past. The molecular analysis is based on N. i. 
caliginosa from western Poland and N. i. in- 
terjecta from the Transcarpathian region of 
Ukraine. 

The situation is intriguing as, in addition to 
the expansion of N. i. caliginosa, also N. i. 
interjecta is expanding northwards within central 
Europe (Nowacki & Pałka 2015, Balint et al. 
2016,). It is therefore highly likely that the ranges 
of these two subspecies will eventually overlap. 
As a consequence, gene pools of the two subspe- 
cies may merge, eventually leading to disappear- 
ance of the two separate taxa. 
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2. Materials and methods 


For the molecular research, adults of both subspe- 
cies of N. interjecta, was accumulated between 
2009 and 2017 during fieldwork carried out at 
Vynohradiv (48°9’3”N; 23°5’8”E) in the Trans- 
carpathian region of Ukraine, and from 2000 to 
2017 in western Poland at the following locali- 
ties: Mrzeżyno (54°8'11”N; 15°1645”E), 
Drawsko Pomorskie (53°31°45”N; 15°48°24”E), 
Unieście (54°16'14°N; 16°5’42”E), Potrze- 
bowice (52°51°5”N; 16°9°53”E) and Poznań 
(52°24’28”N; 16°51°38”E). The moths were at- 
tracted at night to a white screen illuminated by a 
250 W mercury vapour lamp. The individuals 
were picked off the screen and placed in test tubes 
containing ethanol: 


— Noctua interjecta interjecta: 5 exx. Vyno- 
hradiv, Transcarpathia, Ukraine, 3.—5.VIII. 
2016, 

— Noctua interjecta caliginosa: 3 exx. Mrze- 
zyno, Poland, 12.VIII.2016; 1 ex. Drawsko 
Pomorskie, Poland, 10. VIII.2016; 1 ex. Poz- 
nan, Poland, 13.VII.2016. 


For the morphological analysis, a larger sample 
of the moths from both subspecies was trapped at 
Vynohradiv, Transcarpathia, Ukraine (N. i. inter- 
jecta) and Drawsko Pomorskie, Mrzeżyno, 
Unieście, Potrzebowice and Poznań in western 
Poland (N. i. caliginosa). 


2.1. Morphological analysis 


The spans and colouration of both pairs of wings 
were compared in ten individuals from each sub- 
species. The genitalia of five males and five fe- 
males of each subspecies were also studied. The 
genitalia preparations were made according to the 
method described by Nowacki (1995). 


2.2. Laboratory procedures 


DNA extractions from a single leg were per- 
formed using a NucleoSpin Tissue Kit 
(MACHEREY-NAGEL, distribution in Poland 
by AQUA LAB) in line with the standard proto- 
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col. For taxonomic identification a 658 bp region 
near the 5’ end of the mitochondrial cytochrome c 
oxidase I (COI) gene was amplified with the 
primers LEPF1 5’-ATTCAACCAATCATAA- 
AGATATTGG and LEPRI 5’-TAAACTTCT- 
GGATGTCCAAAAAATCA (deWaard et al. 
2008, Hausmann etal. 2011). The amplified frag- 
ment includes the region regarded as the standard 
barcode for the animal kingdom (Hebert et al. 
2003). 

PCR reactions were performed in 30 ul ofre- 
action mix containing 1 x REDTaq® Ready- 
Mix™ PCR Reaction Mix (Sigma-Aldrich, 
Darmstadt, Germany), 0.2 uM of each primer and 
1.5 ul of the extracted DNA. The amplification 
conditions were as follows: initial denaturation 
for 3 min at 94 °C, followed by 40 cycles for 30 s 
at 94 °C, 30 s at 55 °C and 1 min at 72 °C, ending 
in a final extension step of 5 min at 72 °C. The 
Exo-SAP clean-up protocol was used for enzy- 
matic purification of the amplified products. The 
reactions were performed by adding 0.025 ul 
Exonuclease I (20 U/l), 0.25 ul Shrimp Alkaline 
Phosphatase (1 U/l) and 9.725 ul ddH,O to each 
of the PCR products, followed by incubation of 
the samples at 37 °C for 30 min and 95 °C for 5 
min. 

The purified PCR products were sequenced 
on the Applied Biosystems® ABI 3500xL Ge- 
netic Analyzer platform (Applied Biosystems, 
Foster City, USA). The resulting DNA chroma- 
tograms were analysed and aligned using BioEdit 
software v.7.0.5.3 (Hall 1999). Although the de- 
signed primers amplified a fragment of about 658 
base pairs in length, only 600 base pairs were 
clearly readable in all the samples. 


2.3. Statistical analysis 


The number of haplotypes (H), haplotype diver- 
sity (h), nucleotide diversity (7) and mean num- 
ber of nucleotide differences among the 
haplotypes (k) in the overall sample and in partic- 
ular populations were calculated using DNAsp 
5.10 (Librado & Rozas 2009). ARLEQUIN 
v3.5.1.2 (Excoffier & Lischer 2010) was used to 
identify unique haplotypes. The genetic distance 
among haplotypes was calculated using the un- 
corrected proportional distance (p-distance) in 


MEGA version 5 (Tamura et al. 2011). A me- 
dian-joining haplotype network (Bandelt et al. 
1999) was constructed in NETWORK v4.6.1.1 
(Fluxus Technology Ltd.). In addition, a neigh- 
bour-joining tree was inferred. We used a Tamura 
3-parameter substitution model, chosen accord- 
ing to the Bayesian Information Criterion using 
ModelGenerator v0.85 (Keane et al. 2006). The 
COI sequence of N. janthe (Borkhausen, 1792) 
(KX045227, Mutanen et al. 2016) was used as an 
outgroup. Finally, the haplotypes identified in 
this study were compared with homologous se- 
quences of N. interjecta from GenBank (https:// 
www.ncbi.nlm.nih.gov/gen bank). 

Divergence times were estimated using the 
600 bp COI gene in Beast V. 1.7.5 (Drummond et 
al. 2012). The HK Y+G model was implemented, 
and the three codon positions were unlinked in or- 
der to be estimated independently. As there are no 
fossil records of reasonably closely related spe- 
cies that can be used to calibrate the tree topolo- 
gies for this study, we used two conventional mu- 
tation rates for the arthropod mitochondrial COI 
gene taken from the literature: 2.3% (originally 
proposed by Brower (1994), with similar rates 
later independently found in other insects; see be- 
low) and 1.5% (Farrell 2001, Quek et al. 2004) 
per million years (per My). BEAST analyses 
were carried out according to either a strict or a 
lognormal relaxed clock so as to study the effect 
of these two clock models on the resulting esti- 
mates of divergence times. 

According to the strict clock model, we fixed 
a rate of 0.0115 or 0.0075 substitutions/site/mil- 
lion years (corresponding to a divergence rate of 
2.3% or 1.5% per million years, respectively). 
For the lognormal relaxed clock model, we as- 
sumed a normal distribution with a mean of 
0.0115 or 0.0075, and set the standard deviation 
at 10%. 

A Yule speciation process model was selected 
for the tree prior, with a 50% majority-rule con- 
sensus tree obtained from Bayesian analysis of 
the combined molecular and morphological 
dataset as the input topology. Two independent 
MCMC analyses were run for 10,000,000 gene- 
rations, with Markov chains sampled every 1,000 
generations. The analysis results of the two runs 
were combined in LogCombiner v1.7.5 with the 
initial 10% of the trees discarded as burn-in. 
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Fig. 1.—a, b. Noctua interjecta interjecta Hübner, 1803 from Vynohradiv, Transcarpathia, Ukraine, 13.VII.2010 
(a) and Vynohradiv, Transcarpathia, Ukraine, 4.VIII.2016 (b). — c, d. Noctua i. caliginosa (Schawerda, 1919) from 
Unieście, Poland, 2.VIII.2000 (c) and Potrzebowice, Poland, 1.VIII.2008 (d). 


Tracer V. 1.5 (Rambaut & Drummond 2007) was 
used to determine convergence, measure the ef- 
fective sample size of each parameter, and calcu- 
late the mean and 95% highest posterior density 
(HPD) intervals for divergence times. The con- 
sensus tree was compiled with TreeAnnotator 
1.7.5 and displayed in FigTree V. 1.4.0. In addi- 
tion, Bayes factors (Kass & Raftery 1995) were 
calculated using log marginal likelihood with 
stepping stone sampling for a comparison of the 
best fit model between the strict and lognormal 
relaxed clock. 


3. Results 
3.1. Morphology 


The forewings of N. i. interjecta from Trans- 
carpathia in Ukraine have a pale tawny ground 
colour with a faint pattern. The barely visible kid- 
ney-mark and oval are of the same hue as the 
wing ground colour, intermittently edged dark 
brown. The outermost, wavy, cross-line is some- 
what darker than the ground colour, while the 
central cross-lines are broken, faint, here and 


there highlighted by black scales. The hindwings 
are yellow with quite a broad black border fol- 
lowing the outer edge (Fig. 1a, b). The wingspan 
is 36-38 mm. 

Noctua interjecta caliginosa from north- 
western Poland differs substantially in appear- 
ance from the nominative form. The forewing 
ground colour is much darker — dark reddish 
brown, becoming somewhat paler nearer the 
wing base; the pattern is barely visible. The kid- 
ney-mark and the oval are difficult to make out; 
they are of the same hue as the ground colour, and 
only here and there delicately outlined in dark 
brown. The outermost, wavy cross-line is darker, 
almost black, and contrasts with the ground col- 
our. The other cross-lines are faint, sometimes 
picked out in a darker brown, which makes them 
somewhat better visible. The yellow-centred 
hindwings have a very broad black border (much 
broader than in N. i. interjecta) following the 
outer edge of the wing and covering much of the 
wing area (Fig. 1c, d). The wingspan is 32-36 
mm. 

Inspection of the genitalia revealed no struc- 
tural differences between the two subspecies in 
either males or females (Figs 2, 3). 
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Fig. 2. Noctua interjecta male genitalia. — a, b. N. i. interjecta, Vynohradiv, Transcarpathia, Ukraine, 4.VII.2016. 


—c,d. N. i. caliginosa, Mrzeżyno, Poland, 18.VII.2014. 


3.2. Genetic analysis 


A 600 base pair long fragment of the cytochrome 
oxidase subunit 1 (COT) for 10 specimens of N. 
interjecta from Poland and Ukraine was ana- 
lysed: 6 polymorphic sites were found within this 
fragment (4 parsimony-informative sites). All 
variable sites represented transitions. The base 
composition was as follows: A — 16.67%, C — 
46.67%, T — 20.00%, G — 16.67%. 

Three haplotypes were identified among the 
10 sequences obtained. The NImizH1 haplotype 
(GenBank Acc. No: MH172470) proved to be 
clearly dominant in terms of its frequency of oc- 
currence in Poland, being found in four out of five 
individuals from Mrzezyno and Drawsko 
Pomorskie (Poland). Only one haplotype — 
NiImizH3 (GenBank Acc. No: MH172472)—was 
found in five Ukrainian samples from Vyno- 
hradiv (Transcarpathia). The NImizH2 haplotype 
(GenBank Acc. No: MH172471) was found in 
just one sample from Poland (Poznan). There 


were five substitutions separating NImizH1 from 
NImizH3 (p-distance = 0.7%) and two separating 
NiImizH1 from NImizH2 (p-distance = 0.3%) 
(Fig. 4). The mean genetic distance between the 
‘Polish" haplotype group and NImizH3 was 
0.8%. Overall haplotype diversity was 0.644 + 
0.101 with an associated nucleotide diversity of 
0.00437. The average number of nucleotide dif- 
ferences among the haplotypes (k) was 2.622. 
The haplotypes identified in this study were 
compared with homologous sequences of N. 
interjecta from GenBank (https://www.ncbi.nim. 
nih.gov/gen bank) (Table 1). The comparison re- 
vealed seven haplotypes. Two of the haplotypes 
identified in this study (NimizH1 and NImizH3) 
matched sequences from GenBank. Aside from 
Poland, NImizH1 was found only in Germany, 
whereas the Ukrainian NImizH3 also turned up in 
the Netherlands, the U.K. and France (Fig. 4, 
Table 1). The fragment amplified in this study 
was shorter than the sequences deposited in 
GenBank. Hence, some polymorphic sites are 
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missing from the present comparison, namely, 
the sequences HM914226 and KX041259 from 
Germany, indicated by NImizH1 in this compari- 
son; these are two different haplotypes based on 
the alignment of a longer (658 bp) fragment. The 
number of base differences per site (p-distance) 
between sequences did not exceed 1.5%. 
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Fig. 3. Noctua interjecta 
female genitalia. — a. 

N. i. interjecta, Vyno- 
hradiv, Transcarpathia, 
Ukraine, 20.VII.2017. 

— b. Noctua i. caliginosa, 
Poznań, Poland, 
13.VIII.2016. 


The Neighbour-Joining tree (Fig. 5) yielded a 
similar topology to the haplotype network. 
Again, two groups of haplotypes were indicated, 
one comprising haplotypes from Poland and 
NIgbH5, the other containing haplotypes 
grouped around NImizH3. Likewise, an identical 
topology was obtained with Bayesian analysis in 


Table 1. Homologous sequences of Noctua interjecta from GenBank used in comparison with NImiz haplotypes 
identified in this study. Acc. number — GenBank accession number. 








Acc. number Locality Haplotype 
Hausmann et al. 2011: 
HM914226 Germany: Bavaria, Oberbayern, Ruhpolding, Roethelmoos NimizH1 
JF415611 Germany: Bavaria, Upper Bavaria, Mue-Obermenzing, ZSM NIgbH4 
HQ563547 Germany: Bavaria, Oberpfalz, Neustadt a.d. Waldnaab, NIgbH7 
Schwarzenbach, Parkstein-Huetten 
Mutanen et al. 2016: 
KX044477 France: Haute Normandie, Seine Maritime, Doudeville NIgbH4 
KX047693 Netherlands: South Holland, Wassenaar NimizH3 
KX043153 United Kingdom: England, Essex, Caston NiImizH3 
KX043111 United Kingdom: England, Norfolk, Gait Barrows NiImizH3 
KX047097 France: Haute Normandie, Seine Maritime, Varengeville NimizH3 
sur Mer, Bois des communes 
KX041259 Germany: Saxony, Zittau NimizH1 
KX044916 France: Haute Normandie, Eure, Fatouville Grestain NiImizH3 
KX048851 Netherlands: South Holland, Wassenaar NIgbH5 
KX043926 United Kingdom: England, Hertfordshire, Sawbridgeworth, NIgbH6 


Sawbridgeworth Marsh 





ENTOMOL. FENNICA Vol. 30 ° Subspecies of Noctua interjecta in C. Europe 39 


(France Oor = Netherlands ron @- © ae 






NIgbH7 


Fig. 4. Median-joining 
haplotype network of 
Noctua interjecta, 
based on 600 bp se- 
quences identified in 
this study (NImiz) and 
GenBank (NIgb) 
haplotypes. The size of 
a circle is approxi- 
mately proportional to 
the frequency of occur- 
rence of the given 
haplotype in the set of 
22 sequences. The pie 
charts show the propor- 
tions of the haplotypes 
that occur in different 
countries in the studied 
set of 22 sequences. 


NigbH5 


BEAST (tree not shown). The median divergence 
time for NImizH1 and NImizH2 was 0.0950 My 
[95% HPD 0.0049—0.2803] and 0.0388 My [95% 
HPD 0.0012-0.1323] for the 1.5% and 2.3% mu- 
tation rates, respectively. The divergence times 
between the “Polish” haplotypes and the haplo- 
types originating from NImizH3 were 0.6397 My 
[95% HPD 0.2414-1.2015] and 0.4192 [95% 
HPD 0.1571-0.7843] for the 1.5% and 2.3% mu- 
tation rates, respectively. 


4. Discussion 
Morphological analysis of N. interjecta individu- 


als representing populations from two different 
regions of Europe separated by a distance of 


NimizH3 


© mutated position 







NIgbH6 
NigbH4 






NimizH2 


NimizH1 


around 1,000 km showed that N. interjecta 
interjecta from Transcarpathia in Ukraine dif- 
fered distinctly in wingspan and wing ground 
colouration from N. interjecta caliginosa trapped 
in western Poland. Individuals with intermediate 
features were not found in either population of the 
two subspecies. Both male and female genitalia 
of the two subspecies were identical. Hence, there 
should be no morphological barriers to hybridiza- 
tion between these two populations. 

Genetic analysis showed that N. i. interjecta 
from Transcarpathia in Ukraine, where only the 
NImizH3 haplotype was found, differs distinctly 
in this respect from N. i. caliginosa from Poland, 
in which the NImizH1 haplotype was predomi- 
nant. Nevertheless, this is quite a small genetic 
distance, less than 1%. In general, it is recom- 
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NIgbH7 (Germany) 


NimizH3 (Netherlands) 
NImizH3 (UK), 2 exx. 
NimizH3 (France), 2 exx. 
NimizH3 (Ukraine), 5 exx. 


NIgbH4 (France) 
NIgbH4 (Germany) 


NIgbH6 (UK) 


NIgbH5 (Netherlands) 


NimizH2 (Poland) 


NImizH1 (Poland), 4 exx. 
NimizH1 (Germany), 2 exx. 


N. janthe 


Fig. 5. Evolutionary relationships of COI haplotypes (600 bp) of Noctua interjecta estimated using the 
Neighbour-Joining method. The haplotypes identified in this study are preceded by Nimiz, haplotypes 
from GenBank by Nigb. The bootstrap consensus tree was inferred from 1,000 replicates. The percent- 
age of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 repli- 
cates) are shown next to the branches. The tree is drawn to scale, with branch lengths in the same units 
as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances 
were computed using the Tamura 3-parameter method (Tamura 1992) and are given in the units of the 
number of base substitutions per site. The COI sequence of N. janthe (KX045227) was used as an 


outgroup. 


mended that interspecies genetic distance, calcu- 
lated based on COI sequence, should be higher or 
equal than 2%, whereas within-species diversity 
of COI should not exceed 2% (Hebert et al. 2003; 
Ball et al. 2005). Hence, the level of differentia- 
tion, found in this study, clearly reflects within- 
species genetic diversity. Surprisingly, however, 
comparison of our results with the haplotype se- 
quences of N. interjecta from other parts of Eu- 
rope (Hausmann etal. 2011, Mutanen et al. 2016) 
shows that the Ukrainian population has geneti- 
cally far more in common with the populations in 
western Europe (France, the Netherlands, the 
U.K.) than with the Polish populations (Fig. 4). 
The results of this research, together with data 
from the GenBank listed in Table 1 and illustrated 


as a median-joining haplotype network and 
phylogenetic tree in Figs. 4 and 5, indicate that the 
NImizH3 haplotype is widespread and character- 
istic of N. i. interjecta, whereas the NImizH1 
haplotype found in Poland and eastern Germany 
is typical of N. i. caliginosa. These results con- 
firm the hypothesis regarding the phylogeny and 
interrelationships of the two subspecies of N. 
interjecta, which assumes that the division into 
two subspecies took place when the species’ orig- 
inal contiguous range in Europe broke up, giving 
rise to the subsequent long isolation of the two 
populations. The genetic analysis suggests that 
this breakup may have occurred between 640 and 
420 thousand years ago. 

Obviously, the precision of genetic dating 
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(estimating the time of divergence between genes 
or species based on DNA sequences) depends 
closely on molecular clock calibration (Roger & 
Hug 2006). Unfortunately, case-specific parame- 
ters underlying the calibration are rarely avail- 
able, so they have to be replaced by some univer- 
sal estimates. As a result, molecular-based esti- 
mates of divergence times frequently do not 
match expectations from non-molecular data 
(Hedges & Kumar 2004). Such discrepancies 
have raised a number of controversies about di- 
vergence time estimates recovered from molecu- 
lar studies (Hug & Roger 2007). 

The universal mutation rate of COI was also 
used in this analysis: hence, the results only 
roughly reflect the age of divergence between the 
two groups of haplotypes and should be inter- 
preted with caution. Nonetheless, the approxi- 
mate time of divergence seems to correspond 
well with the glaciation history of the continent. 
The range disruption could have come about dur- 
ing the Pleistocene as a result of the San (southern 
Poland) or Elster (Germany) glaciation some 500 
thousand years ago, when the large part of Eu- 
rope, including the whole territory of the contem- 
porary Poland, became submerged under an ice- 
sheet. The remainder of the continent, from the 
Pyrenees to the Black Sea, belonged to the 
periglacial zone with permafrost (Mojski 1993). 

Consequently, the following phylogeogra- 
phic scenario appears likely for N. interjecta. The 
species may have been able to survive in two 
refugia — the Pyrenean one along the Atlantic 
coast, and the near-eastern one to the south of the 
Black Sea. It should thus have been during the 
glacial period that the two subspecies N. i. 
interjecta and N. i. caliginosa diverged. Once the 
ice-sheet had receded from Europe in the Holo- 
cene, both subspecies began to spread north- 
wards. Noctua interjecta caliginosa with the 
NimizH1 haplotype, adapted to the Atlantic cli- 
mate, moved primarily up the coasts of the Atlan- 
tic Ocean and North Sea, thereafter gradually ex- 
panding eastwards to Germany and Poland. 

Noctua interjecta interjecta,;which came into 
being in the near-eastern refuge, very probably 
with the dominant NImizH3 haplotype and 
adapted to the warmer Mediterranean climate, 
was most likely extending its range gradually 
westwards along the coasts of the Mediterranean 


Sea as early as the early Holocene. This subspe- 
cies must have occurred during most of the Holo- 
cene in a contiguous distribution in Europe south 
of the line of central European mountains, i.e. the 
Carpathians, Sudetens and Alps, and also the 
highest mountain ranges in the Balkans. The 
meeting of the two subspecies in western Europe 
led to formation of populations with admixture of 
haplotypes, characteristic of both subspecies, 
which are currently found in the Netherlands and 
south-west Germany. 

At this time, however, N. i. caliginosa with 
the NImizH1 haplotype was already present 
much farther to the east in Germany and Poland. 
The climate in these countries, and also along the 
Atlantic and North Sea coasts, is cooler and too 
severe for the Mediterranean N. i. interjecta, 
which has prevented the latter from colonising 
those regions. In this way we can explain that the 
population of N. i. caliginosa with the dominant 
NImizH1 haplotype occurs in Poland. The 
NImizH2 haplotype, derived from NImizH1 
probably later on, and currently it has not been 
found in western European populations. How- 
ever, analysing the distribution of the subspecies 
in Western Europe is beyond the scope of the 
present paper, and awaits for accumulation of re- 
spective data. 
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